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An Alternative Route to D-wave Superconductivity: the Charge Channel 
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We demonstrate that a d^2_y2 superconducting pairing is naturally associated with the screened 
Coulomb interaction 14// (q,^) within the Cu02 plane of high Tc superconductors. This pairing 
instability arises (via the dielectric constant) from wave- vector structure in 14// (q, Ci)) at q = (tt, tt), 
which is associated with 2D Van Hove and local field effects. Our results, which are independent 
of the bandstructure, suggest that direct Coulomb intereactions are significant (when compared to 
the spin channel) and will, at the least, enhance any other underlying mechanism for d^2_y2 wave 
pairing. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.72.-h, 72.30.-|-q. 



There is a growing body of evidence to support the 
claim that, at least in some of the high Tc cuprates, 
the superconducting order parameter has rf-wave sym- 
metry [1]. This evidence, although by no means conclu- 
sive, is strongest for one particular compound YBCO. 
While, it has been argued theoretically [1] that this d- 
wave state may be driven through exchange of antifer- 
romagnetic spin fluctuations, neutron experiments above 
Tc in YBCO repeatedly fail [2] to observe all but very 
weak antiferromagnetic peaks associated with the opti- 
mally doped material. 

The absence of significant spin fiuctuation effects has 
motivated us to investigate the charge channel as a possi- 
ble source for rf-wave attraction. Other workers have con- 
sidered charge fiuctuation induced superconductivity pri- 
marily via multi-band or excitonic effects, though, within 
this channel, none of these were found to lead to d- wave 
pairing [3]. In this paper we investigate an alternative, 
albeit one band scenario, based on the charge dynamics 
and associated wave vector structure. We demonstrate 
that the screened Coulomb interaction 



^eff{(l,<^) = Vo{q)/e{q,oj) 



(1) 



while repulsive in q space, is attractive in the vicinity of 
nearest neighbor spatial separations in real space. Here 
Vo = 2Tre^ /cgq for a 2D system [4] (and Cg is the static di- 
electric constant of the material). Moreover, this struc- 
ture leads naturally to a dj;2_y2-wave order parameter 
symmetry [5]. 

It was pointed out by Kohn and Luttinger [6] more 
than three decades ago that superconductivity may arise 
from repulsive Coulomb effects. This attraction is as- 
sociated with the Kohn anomaly or singularity in the 
dielectric constant at momentum transfer q = 2kp which 
gives rise to a long range oscillatory interactions (Friedel 
oscillations) in co-ordinate space. Because these studies 
were performed on a three dimensional electron gas, the 
resulting superconducting instability was weak and in a 
high angular momentum channel. 

In order to generalize the Kohn Luttinger picture in 
the presence of a lattice, we follow earlier work by Sham 



[8] and Pandey et al [9] who presented calculations of the 
dielectric constant which included local field corrections. 
The screened Coulomb potential defined in Eq. (1) can 
be written as 

V;//(q+G) = v;(q+G) (2) 

X [l + ^x(q + G,q+G')V;(q + G')] 
G' 

Here G is a reciprocal lattice vector and q is in the first 
Brillouin zone. For notational simplicity we drop the uj 
indices. If we define Xgg' (q) = x(q + G, q -|- G') then 
we have 

XGG' = XgG' + X]^GG""^'o(q+ G")XG"G' (3) 
G" 

within the Random-Phase- Approximation (RPA). The 
noninteracting polarizability is given by 



XGG'(q,^) = 2X1 



Wk - Wk-l-q 
UJ + E\^- -Ek-l-q + 



< k-Fq|e'('i+°') "-|k> 



(4) 



where the basis wavefunctions 



^k(r) =< r|k>= ^ VC7ke"^^^',/-(r-R0 (5) 

are the Bloch wavefunctions and (/)(r — R;) are the as- 
sociated Wannier wavefunctions at lattice site R;. Here 
Ck is a renormalization constant. After some algebra, 
Eq.(4) can be written as: 



o I d^k Wk - Wk-Fq 

XGG'(q,-) - 2 (2.)2c. + i?k-i?k+q + i<5 

X ICkPlCk + qP 



(6) 



where 



1 



REVTeX 3.0 



Levin Group Preprint, 1995 




We apply this approach to the cuprates by incorporating 
standard parameterizations of the band structure in a 
tight binding square lattice: 

= —2t[cos{kx) + cos{ky)] + At' cos{kx) cos{ky) (8) 

where t is the nearest and t' the next-nearest-neighbor 
hopping energy. The Wannier wavefunctions can be fit- 
ted to a Gaussian, whose width is adjusted to corre- 
spond to the known [7] values of the ratio t' /t. One can 
then determine the non-interacting polarizability given in 
Eq.(6), and from this the dielectric constant and effective 
interaction are deduced following Eqs. (2) and (3). 

To place our results in the proper perspective, in Fig- 
ure 1 we plot the wave- vector dependence of the calcu- 
lated, static Coulomb interaction Veff{q) for the stan- 
dard YBCO bandstructure as comparedwith the dynam- 
ical spin susceptibility; the latter is derived from fits to 
low uj neutron experiments [2]. By normalizing both 
curves to the same average value, this figure represents a 
plot of the relative importance of wave- vector structure 
in the spin (dotted line) and charge (solid line) channels 
[10]. Under the assumption that the electrons couple 
with equal strength to each channel, and with the same 
"Debye" frequency, a weak coupling calculation would 
lead to a considerably higher ((ij;2_j,2-wave) Tc for the 
charge induced interactions. 

2.0 \ ' 1 
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FIG. 1. Comparison of the calculated interaction in the 
charge channel (screened Coulomb pairing potential) in 
YBCO (at Ef = —l.6t) (solid line) with that deduced phe- 
nomenologically in the spin-channel (dashed line). 

To understand the important role of the lattice, we first 
compare with the 2D electron gas RPA case, at w = 0, 
where 

e{q,io) = l-Vo{q)x%(l,io) (9) 



and 

X q,^^ =-2/77^ ^ — , p , 10 

J (27r)^ w - E'k+q + Ei^ + iS 

where Eq = h^q^ /2m. In the static limit this latter func- 
tion exhibits a Kohn anomaly as seen in the dashed line 
in the inset of Fig. 2. The resulting screened Coulomb in- 
teraction shown in the main portion of the figure (dashed 
line) starts to rise at q = 2kp, forming a peak at g > 2k p. 
This dashed line was found as a limit of our calcula- 
tional scheme in the case of extremely low electron filling 
(where the electronic dispersion is quadratic). For more 
general values of the electron count, and for the usual 
bandstructure associated with LSCO, the screened inter- 
action, derived from Eqs. (2) and (3) is given by the solid 
line in the main portion of the figure. The solid line 
in the inset of Fig. 2 represents the deduced form for 
the associated x°(q) which essentially includes local field 
corrections, but fitted to the form of Eq.(9). This suscep- 
tibility should be contrasted with the Lindhard function 
which takes the form of Eq.(lO) but in which the inte- 
gration is limited to the first Brillouin zone, so that local 
field terms are absent. It should be noted by comparing 
Fig. 1 and Fig. 2 that there only small differences be- 
tween the behavior of the screened Coulomb interaction 
in the two different cuprates. 
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FIG. 2. Comparison of the effective Coulomb potential in 
LSCO at the Van-Hove point (solid line) with the case of a 
2D free electron gas ( dashed line). The upper inset plots 
the effective susceptibility for these two cases, illustrating the 
effect of the Van-Hove singularity on the Kohn anomaly. The 
lower inset demonstrates how Van-Hove effects are enhanced 
by Umklapp processes. 

The principal difference between the full lattice calcu- 
lation (solid line) and the electron gas limit (dashed line) 
can be seen from Fig. 2 to be associated with both the 
small q and q = (tt, tt) regimes. The latter is the more 
important, since the former is associated with very long 
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distance Friedel oscillations, which are not expected to 
play a role in the superconductivity. The origin of the 
structure around (tt, tt) in both Fig. 1 and Fig. 2 is a 
"two particle" Van- Hove effect which may be viewed as 
a "depletion" in the intensity of the screened interaction 
(relative to the dashed line) as a consequence of a max- 
imum in the diectric function. This in turn arises from 
the fact that in a 2D lattice, the density of states in mo- 
mentum space has divergences at the Van-Hove points of 
the Brillouin zone edge. For a square lattice (appropriate 
to the cuprates), these points are at (±7r, 0) and (0, ±7r). 
In this way when the momentum differences between the 
initial and final states in an electric field induced po- 
larization process correspond to the separation of two 
Van-Hove points, one expects a peak in the dielectric 
constant. Since most Fermi surfaces do not necessarily 
intersect the Van-Hove points, the peak position is de- 
termined by the momentum transfer between two points 
on the Fermi surface which are in the neighborhood of 
the Van-Hove points. Here we see that local field correc- 
tions which enter via Umklapp processes are important 
so that the "two particle" Van Hove effect necessarily in- 
cludes contributions from momentum space beyond the 
first Brillouin zone. These processes, which enhance the 
Van Hove contributions further, are shown in the lower 
inset of Fig. 2 [11]. 

The above results lead to important implications for 
the Kohn Luttinger model for superconductivity, when 
generalized to a 2D lattice. It should be stressed that 
the original calculations were based on a weak coupling 
model, as well as the assumption that Migdal's theorem 
is valid for electronic mechanisms. The latter is a highly 
complex issue which we will not adress here [12] since the 
issue of primary interest is the order parameter symme- 
try rather than magnitude of Tg. For these approxima- 
tions it is relatively straightforward to deduce the sym- 
metry of the order parameter. In the case of an electron 
gas, the screened Coulomb interaction depends only on 
|q|. Because there is no directionality, there can be no 
(ij;2_j,2-wave instability. By contrast, in the 2D lattice, 
the screened interaction depends on the wave-vector di- 
rection. Moreover, there is a maximum slightly above 
(tt, tt) as shown in the first two figures. In this way we 
consistently find that the resulting superconducting or- 
der parameter has dj;2_y2-wave symmetry independent of 
the electronic filling and the band structure . A typical 
gap function is shown by the solid line in the inset of Fig. 
3. 

An issue of primary interest in understanding the pair- 
ing mechanism is determining an associated dynamical 
energy scale, which can lead to sufficiently high Tc- This 
is all the more problematic in the context of the low 
values of the boson coupling constants which appear in 
transport experiments [13]. For the spin fluctuation case, 
the energy scale observed from neutron data [2] is com- 
parable to that of phonons, and thus it may not be suf- 



flciently high. 
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FIG. 3. Inverse of dielectric constant for YBCO (at 
Ef = — 1.6t) as a function of frequency: solid line lo = 0; 
dashed line, lo = 3.2t; dash-dotted line, lo = 6.4t. The corre- 
sponding superconducting gaps of the solid and dashed cases 
are shown in the lower inset. Upper inset plots the dispersion 
relation for the plasmon mode (solid line) and the edge of the 
particIe-hoIe continuum (dashed line). 

At flnite frequency, there are two types of relevant ex- 
citations contained in the screened Coulomb interaction. 
One is the particle-hole "exciton" or continuum and the 
other, the usual plasmon collective mode (uip) . The di- 
electric constant vanishes at the latter while its inverse 
passes through zero at the continuum edge, (wea;)- 
As shown in the upper inset of Fig. 3, there are two 
regimes in the parameter space called I, II. At flxed q, 
the flrst satisfles ujp > uigx and the second ujp < uigx- The 
dashed line of Fig. 3 shows the behavior of for mod- 
erate, real frequencies ( uj = 3.2t) which correspond to re- 
gion I. While there is a sign change at small wave- vectors, 
the feature at (tt, tt) still dominates the interaction and 
the resulting (flnite frequency) order parameter is essen- 
tially the same as for the static limit, as can be seen from 
the lower inset. At higher frequency (uj = 6.4t) , corre- 
sponding to II, the behavior is plotted as the dot-dashed 
line of Fig. 3. For these high frequencies, the d wave 
state is no longer stable. In this way one deduces a char- 
acteristic pairing energy for d wave superconductivity, 
which corresponds to that frequency where the plasmon 
mode enters the particle-hole continuum. For YBCO, we 
estimate this cutoff frequency uJcuto} / ~ 4t, which for rea- 
sonable values oft ranges from 250 to 1000 meV. Clearly 
a full dynamical Eliashberg calculation should ultimately 
be performed to incorporate this flnite frequency effect, 
although this is also subject to earlier concerns about 
Migdal's theorem. Nevertheless the present calculations 
demonstrate that the energy scale for the charge channel 
should be considered to be of the order of hundreds of 
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meV, as contrasted with tens of meV in the spin chan- 
nel. Furthermore, this behavior appears consistent with 
Uemura's scaling plots [14]. 

Finally, we note that the presence of q structure in 
e(q,w) also has important implications for the screened 
electron-phonon interaction and its contribution to any 
superconducting pairing. Assuming that the phonon 
propagator is relatively isotropic so that it can be ap- 
proximated by a coupling constant of strength Xp, the 
total effective interaction is then 

Veff = Vo/e-Xp/e (11) 

While, in principle, pseudopotential effects should be 
included, we may incorporate these phenomenologically 
through an effectively stronger Xp. In Fig. 4, we present 
the calculated superconducting gap for various values of 
the phonon coupling strength. The figure illustrates the 
extreme case where the Fermi energy is close to the Van 
Hove singularity, so as to emphacize Van Hove effects. In 
the absence of phonons, as well as for the case of small 
Xp = 0.5 , a similar dj;2_y2-wave gap is obtained. With 
increasing coupling the gap first evolves to a twelve node 
d-wave state { Xp = 1), then to a highly anisotropic s- 
wave state ( Ap = 2) and ultimately ( Ap = 5) to a some- 
what more isotropic s-wave symmetry. The Van Hove 
singularity plays an important role in, whenever possi- 
ble, maximizing the gap values at the four Van Hove 
points. Somewhat away from the Van Hove singularity, 
we find that an intermediate d^y state arises before the 
transition to s-wave symmetry. The behavior discussed 
here should be contrasted with the complementary cal- 
culations in Ref. [15], where the role of Van Hove effects 
on the electron phonon interaction were also studied [16]. 
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FIG. 4. Superconducting gap for YBCO at the van-Hove 
point for various values of the electron- phonon coupling 
strength (Xp). 

In summary, we find important wave vector structure 
in the dielectric constant and therefore in the screened 



Coulomb interaction at q = (7r,7r). This structure is 
associated with the 2D Van Hove singularities in con- 
junction with local field corrections, and leads to dj;2_y2 
wave pairing. This physical picture should be viewed as a 
charge analogue of the standard scenario used to deduce 
d-wave superconductivity from the spin channel. It is the 
amplitude and position of q structure in the interaction 
which determines the nature and degree of attraction in 
real space. Future calculations will extend this work to 
the multi-band case. Whether this provides all or only 
a portion of the superconducting attraction, our calcula- 
tions demonstrate that the direct Coulomb channel will 
always act to enhance any (other) underlying propensity 
for d- wave superconductivity. 
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